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Abstract 

We investigate zero-field Ising models on periodic approximants of planar quasiperiodic tilings by means of partition 
function zeros and high-temperature expansions. These are obtained by employing a determinant expression for the 
partition function. The partition function zeros in the complex temperature plane yield precise estimates of the critical 
temperature of the quasiperiodic model. Concerning the critical behaviour, our results are compatible with Onsager 
universality, in agreement with the Harris-Luck criterion based on scaling arguments. 
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1. Introduction 

The recent, rather controversial discussions about 
magnetic ordering in quasicrystals, in particular for 
alloys containing rare earth elements jl], |^, has re- 
vived the interest in simple models of magnetic or- 
dering in quasicrystals and their properties. The 
influence of an aperiodic order on magnetic phase 
transitions can be understood by heuristic scaling 
arguments Q resulting in the Harris-Luck cri- 
terion for the relevance of aperiodicity. For one- 
dimensional aperiodic Ising quantum chains, equiv- 
alent to two-dimensional layered Ising models with 
one-dimensional aperiodicity, exact real-space renor- 
malization techniques can be used that prove the 
heuristic criterion for a large class of models based 
on substitution systems [|[ . For the classical Ising 
model defined on two-dimensional quasiperiodic cut- 
and-project tilings, the criterion predicts that the 
critical behaviour is the same as for periodic lattices, 
hence these models belong to the Onsager universal- 
ity class. This has unanimously been corroborated 
by results obtained from various approaches |7| , and 
has been commonly accepted. 

Nevertheless, it is interesting to study these mod- 
els in more detail, and in this short note we present 
results for periodic approximants for the Penrose and 
the Ammann-Beenker tiling. The zero-field Ising 
model on periodic approximants can be treated by 
methods that allow us to compute expansion coef- 
ficients and partition function zeros for the infinite 
periodic tiling, hence we do not have to deal with 



finite-size corrections as one usually accounts if only 
finite patches are considered. Instead, we can study 
the dependence on the size of the unit cell. From the 
partition function zeros, we obtain, in principle, the 
exact critical temperatures for the periodic approxi- 
mants. In this way, we are able to derive precise esti- 
mates for the critical temperatures of the quasiperi- 
odic models. For the high-temperature series of the 
free energy, we investigate the fluctuations in the se- 
quence of expansion coefficients observed in |^ , using 
the exact values for the critical temperatures of the 
periodic approximants. 

2. The Ising model on periodic approximants 

We place an "Ising spin" ai — ±.l on each vertex 
i of the periodic approximant. A configuration <t = 
{a I, (72, cr„) is assigned the energy 
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where J > is the ferromagnetic exchange coupling, 
which we choose identical for all pairs of neighbour- 
ing spins, and we sum over all neighbouring pairs 
of spins {i,j) at vertices i and j. Two vertices are 
neighbours if they are connected by an edge of the 
tiling. 

There exist several variants for the solution of the 
zero-field Ising model; the method we use here was 
initially formulated by Kac and Ward for the square 
lattice Ising model M. The Kac- Ward determinant 
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Fig. 1. Part of the partition function zeros in the complex 
z plane for the third periodic approximant of the Penrose 
tiling. 

and its application to periodic approximants was de- 
scribed in detail in our recent paper Let us just 
mention briefly that the partition function Z{f^) — 
exp[— /3i5(cr)] for a zero-field Ising model on an 
arbitrary planar graph, where /3= X/ksT denotes the 
inverse temperature, is expressed as the square root 
of a determinant of an infinite matrix K(w), where 
w = tanh{pj) For a periodic approximant, this 
matrix is cyclic. Thus the problem can, by Fourier 
transform, be reduced to the computation of a deter- 
minant of a finite matrix K{w, (pi, (^2) that depends 
on two additional parameters <^i, V'2 G [0,27r) which 
play the role of Bloch phases. The size of the ma- 
trix K(w, ipi, (P2) is given by the number of oriented 
edges in the unit cell, hence by AN for the Penrose 
and the Ammann-Beenker approximants, where N 
denotes the number of spins in the unit cell. Apart 
from a factor —1/(5, the free energy per vertex in the 
thermodynamic limit /(w) —\miM~too\^o^Z{P)/M] 
is given as 



where M denotes the the number of spins in a finite 
patch approximating the infinite system. 

For large approximants, it is expendable, and not 
very profitable, to compute the determinant analyt- 
ically as a function of the parameters w, ifi, and ip2- 
Instead, we rather calculate the partition function 
zeros or the coefficients of the high-temperature ex- 
pansion of the free energy by exploiting the fact that 
the matrix K has the form 

K{w,(pi,ip2) = I + wL{ipi,ip2) (3) 



Fig. 2. Same as Fig. but for the second periodic ap- 
proximant of the Ammann-Beenker tiling. 

where I denotes the unit matrix. In particular, the 
matrix L does not depend on the variable w. This 
makes it possible to compute the partition function 
zeros for given values of (pi and Lp2 from the eigen- 
values of the matrix L(tpi,ip2), and the complete set 
of partition function zeros is just the accumulation 
of the zeros for all ipi,(p2& [0, 27r). The coefficients 
of the high-temperature expansion of the free energy 

00 

fM = E52„^'" (4) 
can be calculated from traces of powers of L by [|j 



which follows by expanding the logarithm in Eq. (|^) . 

3. Partition function zeros 

In what follows, we adopt the usual convention and 
represent the partition function zeros in the variable 
z = {l + 'w)/(l — w) = exp(2/3J). For an infinite sys- 
tem, the zeros of the partition function in the com- 
plex z plane accumulate on lines or areas that sepa- 
rate different regions of analyticity of the free energy, 
hence different phases of the model. In particular, 
partition function zeros on the real axis correspond 
to phase transition points, and thus determine the 
critical temperature. The behaviour of zeros around 
these points is governed by the corresponding critical 
exponents |Tl| . 

From the present approach, it is obvious that the 
zeros themselves are parametrized by the two angles 
Lpi and Lp2, thus we expect that they generically fill 



Table 1. Critical temperatures Wc = tanh ( J/fcsTc) for 
periodic approximants of the Penrose tiling, extrapolated 
to the quasiperiodic case with an estimated error. Here, 
m labels the approximants with A'' spins per unit cell. 



m 


AT 


Wc 


1 


11 


0.401440 380 


2 


29 


0.395 411099 


3 


76 


0.395 082 894 


4 


199 


0.394 554 945 


5 


521 


0.394 523 576 


6 


1364 


0.394454 880 


7 


3571 


0.394451035 


8 


9349 


0.394441450 


9 


24476 


0.394439 826 


10 


64079 


0.394439 319 


oo 


oo 


0.394 439(1) 



areas in the complex plane. In the square lattice 
case, however, the zeros fall on two circles with ra- 
dius v^, centred at z = ±l; similar results are found 
for other simple planar lattices. In order to illus- 
trate the behaviour for the periodic approximants, 
we show the zero patterns for the third approximant 
of the Penrose tiling with N = 76 spins per unit 
cell in Fig. ^, and the result for the second approx- 
imant of the Ammann-Beenker tiling with = 41 
in Fig. 1^. Here, the zeros were computed for angles 
(pi,(^2e{m7r/20 I TO = 0,1, 2,..., 39}. 

Clearly, the patterns are more complicated than 
the two circles found for the square lattice, in partic- 
ular further away from the positive real axis. Wc find 
two zeros on the positive real axis, one for ipi=ip2 — 
and the other for ipi = ip2=TT, corresponding to the 
ferromagnetic and antiferromagnetic critical points, 
respectively, which have the same properties due to 
the bipartiteness of the tilings. In the vicinity of the 
critical points, the zero patterns are very well de- 
scribed by segments of circles that orthogonally in- 
tersect the real axis, showing that the corresponding 
critical exponent a stays the same as for the square 
lattice, i.e. a = 0, corresponding to a logarithmic sin- 
gularity of the specific heat. The zeros on the real 
axis determine the critical temperatures for the peri- 
odic approximants, the results for the ferromagnetic 
critical point are shown in Tables |l| and |[ The values 
are in agreement with results of recent Monte Carlo 
simulations using the invaded cluster algorithm [|l2j . 
Interestingly, the values for the Ammann-Beenker 
approximants in Table appear to converge much 
faster than those for the Penrose case in Table 
This might be related to the observation that the 
mean coordination mmrbers jl^ , in particular the 
number of next-nearest neighbours, also converge 



Table 2. Same as Table y, but for periodic approximants 
of the Ammann-Beenker tiling. 



TO 


N 


Wc 


1 


7 


0.396 850 570 


2 


41 


0.396 003 524 


3 


239 


0.395 985 346 


4 


1393 


0.395 984 811 


5 


8119 


0.395 984 795 


oo 


oo 


0.395 984 79(2) 



faster for the Ammann-Beenker case, thus the pe- 
riodic approximants of the Ammann-Beenker tiling 
are, in this sense, "closer" to the infinite quasiperi- 
odic case than those of the Penrose tiling. 

4. High-temperature expansion 

In principle, the expansion coefficients of the high- 
temperature series can be calculated exactly for 
quasiperiodic cut-and-project tilings in the frame- 
work of the projection method. However, the num- 
ber of graphs contributing to the expansion grows 
tremendously, thus limiting the applicability of this 
method. Recently j^, we presented the leading 
terms, up to 18th order in the expansion variable 
w = tanh(/?J), of the high-tcmpcrature series 
for the free energy f{w) of the Ising model on the 
Penrose and the Ammann-Beenker tiling. The se- 
ries coefficients contain information about the criti- 
cal behaviour of the model. Assuming that the free 
energy f{w) at the critical point Wc shows a power- 
law singularity f{w) ~ (l — w'^/wl) , the ratios of 
successive coefficients behave as 

-^-^(l-^) (6) 

for large n, compare Eq. (^). Thus, plotting the ra- 
tio g2n/ 92/1-2 as a function of 1/n, the data points 
should, for sufficiently large values of n, lie on a 
straight line. From its intercept at 1/n = and 
its slope, we can, in principle, extract the values of 
the critical temperature Wc and the critical exponent 
K = 2 — a ~ vd. Here, d — 2 denotes the spatial di- 
mension, V is the critical exponent of the correlation 
length, which is = 1 for the Onsager universality 
class, and the last equality is due to scaling. 

As shown in this plot reveals large fluctuations, 
in particular for the smallest periodic approximants, 
while it behaves very nicely for simple lattices as, 
for instance, the square lattice. These fluctuations 
rendered impossible any reasonable estimate of the 
critical temperature or the critical exponent from the 
series expansion, at least in this way. In order to get 
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Fig. 3. The differences r^n (|7|) between the ratios of ex- 
pansion coefficients and their expected asymptotic be- 
haviour (|^ for the three smaiiest periodic approximants 
of the Penrose fifing. For cfarity, the data for m = 1 and 
m = 2 were shifted. Lines are incfuded to guide the eye. 



some more insight into tfie nature and tfie reasons of 
the fluctuations, we extended the series expansions 
for the periodic approximants considerably. 
In Fig. 1^, we plot the differences 

= (7) 



n wf. 



92n-2 



between the ratios of the expansion coefficients, up 
to order 2?!= 128, and their expected asymptotic be- 
haviour (j^) , with K = 2, for the smallest approximants 
of the Penrose tiling. Here, we use the critical tem- 
peratures Wc given in Table ^ Apparently, the ratios 
r2n for the smallest approximant with A'^^ll spins 
per unit cell fluctuate strongly for 20 ~ 2n ~ 60. For 
larger values of 2n, the deviations from the asymp- 
totic behaviour are small, similar to the square lat- 
tice case. The second approximant with = 29 
shows a strikingly similar region with strong fluctu- 
ations, but at larger orders, starting around 2?! ~ 50 
and waning just around the maximal order 2n = 128 
in Fig. ^. For the third approximant, no comparable 
region of particularly strong fluctuations is observed, 
but we might speculate that the order where these 
fluctuations occur grows roughly with the number of 
spins in the unit cell. In this case one has to go to 
higher orders to observe the effect. One might expect 
that for orders clearly beyond the number of oriented 
edges in the unit cell, i.e., for n^4iV, the deviations 
from the asymptotic behaviour are small. 

For larger approximants, the expansion coeffi- 
cients g2n rapidly approach those of the quasiperi- 
odic tilings Thus, we may draw the conclusion 
that the data for the periodic approximants show 
cross-over phenomena between the characteristic be- 
haviour for the quasiperiodic model and the simple 
square-lattice case, although for the cases considered 
here the asymptotic behaviours should coincide. It 
seems reasonable that the cross-over occurs approx- 



imately at the size of the unit cell. 

5. Concluding remarks 

The Ising model defined on quasiperiodic tilings 
such as the Ammann-Beenker and the Penrose tiling 
belongs, according to scaling arguments, to the same 
universality class as the model on the square lattice. 
This is in accordance with our results on the partition 
function zeros and the high-temperature expansions 
for the corresponding periodic approximants. Nev- 
ertheless, the quasiperiodic systems show interesting 
features, and a better understanding of the cross- 
over phenomena in the high-temperature expansions 
of the periodic approximants is desirable. 
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